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INCLUSION OF KNOWN INTEGRALS IN THE OPTIMAL 
TRAJECTORY PROBLEM 

INTRODUCTION 


The optimal trajectory problem concerning a single attracting center — which is 
usually taken to be that trajectory which minimizes mass expenditure — has been treated 
in numerous publications, 

Two important integrals of the problem were quickly isolated in the early 
literature. One of these integrals, vectorial in nature, corresponds to a sort ol “angular 
momentum” conservation. The other, a scalar integral, is a sort of conservation ot 
“energy.” It has long been realized that, if full use could be made- of the information 
contained in these integrals, a great reduction in the number of numerical integrations of 
the adjoint equations would be achieved. Thus, the iteration required would be faster. 

The transformation which yields such a reduction has been elusive. In this paper 
the transformation is explicitly demonstrated. The transformation exacts a price, however, 
in that the resultant equations are ill-behaved from the numerical point of view. 

In the case of a planar trajectory which has the magnitude of the radius and 
velocity vector, as well as the flight path angle, constrained at the terminal point, a 
regularization is achieved. In the general three-dimensional case, such a regularization has 
not been found. 

The first section below deals with known results. Subsequently, the planar case is 
discussed in depth and the formulation of the three-dimensional case follows. A few 
concluding remarks on other unsuccessful lines of approach to achieve a general 
regularization are given. 


KNOWN RESULTS 


The general equation for rocket motion about a single attracting center is 



( 1 ) 


where "r is the position vector of the rocket from the attracting center, T is the vectorial 
thrust, m is the mass of the rocket, and ju is the gravitational parameter of the attracting 
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center. Derivatives with respect to time are shown by the dot notation. The rate of 


change of the mass is related to the magnitude ot the thrust, T, and ihe exhaust velocity, 
c, via 


m = - T/c 


( 2 ) 


If one defines X as the solution of the vectorial differential equation 



(3) 


and further defines 


S 



a 

c 


(4) 


where 



(5) 


then it is well known [1] that an optimal trajectory (in the sense of minimizing mass 
expenditure) will result if the direction and magnitude of thrust are chosen according to 


T 



( 6 ) 


^max ^ S > 0 


T = 


T mi„ if S < 0 


(7) 


The case where S = 0 for a nonzero interval of time will not be considered here (see 

Reference 2). 
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Two important integrals of the system of equations represented by equations (1) 
through (6) are given by 


r - X 


r = C 


( 8 ) 


and 


TS - i = c 4 

T 3 


(9) 


where C = (C, , C 2 , C 3 ) and C 4 are integration constants. 

If l?|, fr|, and 6 ^0 = cos' 1 are fixed at the final point and only these 
quantities are fixed, then it may be shown that f 1 ] 


C = (C 15 C 2 ,C 3 ) = O = (0,0,0) 


( 10 ) 


Furthermore, if the final time is unspecified, it is known that [11 


C 4 = 0 


(ID 


THE PLANAR CASE 


The case represented by equations (10) and (11) will be discussed first. Assuming 
that equation (10) is valid,- take the dot product of r with equation (8) and interchange 
the order to yield 


X • (7 X t) = 0 . (12) 

Geometrically, this result can be interpreted to indicate that the X vector will 
always remain orthogonal to a vector which is orthogonal to the plane formed by ~r and "r! 
Thus, one can write 
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( 13 ) 


£ = 


-» 

£ r + r? v 


where £ and tj are scalar functions of time. The degrees of freedom represented 
by £ and t? are utilized by requiring that equations (8) and (9) be satisfied identically. 
In Reference 1 it is shown that 


£ - 




(14) 


and 


V 



- 2 £ 


(15) 


where 


-*• 

v 



(16) 



(17) 


The numerical application of the theory developed to this point is now clear. 
From appropriate starting values of the required quantities, equations (5), (14), and (15) 
could be integrated forward in time. Equation (13) then yields the X vector as a 

function of r and v, so that equations (6) and (7) give the direction and magnitude of 
the thrust vector. Equation (1) may then be integrated, numerically, to yield updated 

— y _ y 

values of r and v . 

The advantage of using equations (14) and (15) in lieu of integrating equation (3) 
directly is that of dimensionality. Specifically, equation (3) represents three second-order 
equations which are not necessarily simple to integrate. Equations (14) and (15) are two 
first-order equations — i.e., a reduction from 7 to 3 in the number of integrations necessary 
to determine the steering program has been achieved. Furthermore, the transversality 
conditions [equation (8)] are identically satisfied; likewise, the fact that the Hamiltonian 
[equation (9)] will remain identically zero is assured. 
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The advantage of reducing the number of required integrations becomes even 
more important in light of the two-point boundary value problem associated with the 
trajectory optimization procedure. The satisfaction of the four equations represented by 
equations (8) and (9) is guaranteed no matter what initial values are chosen tor the 
variables £ and r? . Thus, a vast reduction in iteration accompanies the reduction in 
integration. 

The difficulty in using equations (14) and (15) is^that the occurrence of r • v in 
the denominator of equation (14) is a near-singularity, r • v will become sufficiently small 
that overflow is almost certain to be encountered during the optimal trajectory 
calculations. Therefore, the next process is the regularization of equation (14). 


Set 


£ - <t>n 


(18) 


Then 

X = y/cj) 2 r 2 + v 2 + 2(p r • v 

Define, for convenience, 

X* = y/<l> 2 r 2 + v 2 + 24> r ■ v 
Using equations (18) and (20) in equation (15) gives 



Next, differentiate equation (18) to yield 

- /• , T0 \ 

£'=0t? + 0t 7=U-2</> 2 - — Jh 
where equation (21) was used for r? . 


(19) 


( 20 ) 


( 21 ) 


( 22 ) 
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Substituting equation (22) into equation (14) and employing equation (18) to 
eliminate £ gives 

»-£ )-ix 

r / m 

The expression on the right side of the first equality sign comes from separating S into 
its elements via equation (4). Subsequently, equation (2) was employed. 

To simplify the algebraic calculations, abbreviate 

$ = (?-v)(0 - 2<p 2 ) - (p(v 2 - 7) - (24) 

so that equation (23) can be written as 


( 

— 


1? (r • v) (0 - 2 <p 2 ) - <p (v 



17 = rh 0 


(25) 


Since the singular arc case has been excluded via choice, the thrust magnitude 
must be either T max or T min . In either case m is a constant. Thus, equation (25) can 

be differentiated to give 


r) <t> + 77 <t> 


m 0 = m- 


TX*t? 


nr 


(26) 


from equation (5). 

Eliminating 17 by use of equation (21) and assuming that 


17 0 


(27) 


allows equation (26) to be written as 





<t> = rh- 


T \* 


m 


(28) 
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The next task is to calculate 4> from equation (24). This calculation is somewhat 
complicated but the result is 


4 > = 


( w 2 _ H\ + _L. (0 r 2 + ~r v )1 (0 - 2 0 2 ) + (r’-'v) (0 - 4 0 0) 

_\ r / m X* J 




-»• -»• 

T ^ 2 , T ‘ v 

2-— (0r-v + v 2 ) - 
m X* r 


T m X* 


mX* 


0 0 r 2 + 0 (?• v) + 0 2 (r • v) + 0 


2/?.^ + A^v 2 - 


m‘ 


T -»■ IX r • v 

+ — (0 r + v) 2 -5— 


m X* 


(29) 


Inserting equations (29) and (24) into equation (28), gathering, and simplifying yields 


(r’.'v) 0 - 0 (?• v)^6 0 + + 4 0 3 (?•■?) - 0 2 ^3?- v + 2 0 r 2 ^ 


mX* 


+ 01 


2TX* 2Tv 2 + M i_v \ _L f £ 141 ) = 0 


m 


m X s 


m X* \ r 


(30) 


Further simplifications of equation (30) are possible. Notice that 


2 T 2 T 2 _ 2T 

X* — v 2 


m 


m X* m X s * 


0 (0 r 2 + 2 r • v) 


(31) 


Then equation (30) may be written as 


0 - 0 l6 0 + 


T \ , T M0 ^ T n 

) + 4 0 3 + — — 0 2 + — r- + — — - 

J m X* r 3 m X* r 


m X* 


(r - v) = 0 


(32) 


7 



In general r • v =£ 0 so the first factor of equation (32) must vanish. Then 
equation (32) becomes 


0-500 + 3 0 3 + f 0 + f - 0 + 0 2 + -p j - 0 . (33) 


This is the equation which was originally sought. There is no longer a singularity 
involving r • v . 

The substitution represented by equation (18) also converts the state equation (1) 
into a very convenient form. Using equations (18), (19), and (20) in equation (1) yields 


r “ 


T 

m X* 


ju r 

(0 r + v) j 

r 


(34) 


Equations (33), (34), and (5) now constitute a complete set of differential equations for 
the optimal trajectory problem. 

Equation (34) may be put into an interesting form as follows: 


r x r - — — v x r , (35) 

m A* 


but 


r x r 


d -* -* 


(36) 


Writing 


L = 


r x v 


(37) 
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equation (35) becomes 


->■ 

L = 


m X* 


( 38 ) 


This equation easily integrates to yield 


L = 


/ 


L 0 e 



(39) 


where L 0 is the initial value of the angular momentum. Thus, L is always parallel to L 0 
(which was known anyway) and equation (39) can be written in scalar form as 


f 


L = L 0 e 


in X : 


dt 


(40) 


A modification of equation (33) which makes the problem more numerically 
tractable follows [31 : 

Set 


0 = P 


(41) 


Then, from equation (33), 


p = 0= : 50p-3 0 3 + 


(* + ( p ‘ ' ?) 


(42) 


Thus 
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or 


d 0 P 




( 44 ) 


Next, set 


0 = 


(45) 


and 


P = 


w e 


2v 


(46) 


Then 


dp 

d0 


d w 
d v 


+ 2we l 


5 e" - 3 


*3 v 


w e 


2v 




m X 4 



e 2v 


w e 


2v 



(47) 


Clearing and collecting yields 


d w 

w + 2 

d 


( w -f) <»-•>-(' + ^ x* - * - *v) 


(48) 


For the trivial case of field-free coast, equation (48) solves to yield 


v - In C 5 


(w - 1 ) 1/z 
(w - 2 ) 1/4 


( 49 ) 


where C 5 is an integration constant. 
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THE NONPLANAR CASE 


Although the initial set of equations, equations (1) through (9), apply equally 
well to either two- or three-dimensional trajectories, equation (12) involves the 
assumption of planarity. The next goal is to relax this assumption by choosing a nonzero 
value for Even so, the time-open problem (C 4 = 0) will be assumed. 

In order to allow for nonplanar steering, write 


X = £r > +7jv'+f(r’ xv*) 

where £, p, and f are scalar functions of time. The variables £ and rj as used in this 
section agree with the variables £ and tj ot the last section if f = 0. From equation (50), 

X* = £ r* + £"v + + 7? ~v + f (T x ~v) + f (r x v) (51) 


and, from equations (1) and (50), 



+ 17 v + f (r x v) 



(52) 


Calculating "rx "v and inserting the result along with equation (52) into equation (51) 
yields 



m + U1? 

r 3 m A 



V + t + “~T 172 

m A. 



+ 




(r 


X 


V) 


( 53 ) 


It is now possible to calculate X x r as 



while from equation (50) 


X x v = £( r x v) + f[(r • v) v - v 2 r] 


(55) 


Using equations (54) and (55) in equation (8) now gives 

+ + • ] 4rr) <7xt) + [( J + ^’ ,f ) r ’ 

f + 7 ] n (7 • v) + v 2 

m X / 

If the dot product equation (56) with 7x 7 is calculated, one finds 




V + 2£ + — (t? 2 - f 2 r 2 ) = 
m X 


C ■ (v x 7) 
l7 X 7 P 


(57) 


Dotting equation (56) with v yields 


t 


+ 


2T 



c ■ 7 
i7 x 7i 2 


(58) 


Finally, dotting equation (56) with r yields 
C • 7 

r = _ _ 

lr x v I 2 


(59) 
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Since two expressions which govern the variable f have been found [equations 
(58) and (59) J , it is necessary to show that they are compatible. This is most easily 
accomplished as follows. From equation (52) 


— ^ d — ^ ^ T — v — j. — ^ — k 

rxv = — (r x v) = — - [tj (r x v) + ? r x (r x v)] 
dt m X 


(60) 


so that 


/-> d -»• -> T -+ -»■ 2 

(r x v) • — (r x v) = — -tj (r x \) 1 
dt m X 


(61) 


Then calculate 


_d 

dt 


1 


(r x v) 2 


(r x v)'* 


(r x v) • — (r x v) 
dt 


- 2 n t 

m X (r x "v) 2 


(62) 


Assuming equation (59), one finds 


C ■ v -> ->d 
t = — — + C • r — 

(?x v) 2 dt 


1 


(r x v) 2 


C • v 2 rj TC • r 


(r x v) 2 m X (r x v) 2 


(63) 


Substituting for the combination of variables represented by equation (59) gives 


2T 
rn X 


r?r 


C • V 
(T x V ) 2 


Thus equations (58) and (59) are compatible. . From this point onward the discussion 
deals only with equation (59) since equation (58) is redundant. 

It is still necessary to derive the equation governing £ . To do this, return to 
equation. (9). Calculating T-v from equation (53) yields 
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/ T r/n Tf 2 ^ ^ ^ /. T 

X • v = ( £ + — - £ 77 - ^ + — - r • v j (r • v) + j 77 + £ + — - v 2 

\ m X r 3 m X / \ m X 


T t 2 

j 

m X 



(64) 


From equation (50), 


X • r = £ r 2 + 17 r 


(65) 


Then equation (9) may be written as 


[i + -x 

L rn x 


T -> -> , C - (r x v) 

. £ 7 ? + — - f 2 (r • v) (r • v) - £ v 2 + _ 

X m X J |7 x 7l 2 


+ £- = TS 


( 66 ) 


• -t 

where equation (57) was used to eliminate r? . Assuming that r • v ^ 0 one finds 


£ - £ 


/b - a\ 

r 


C ■ ( r x v) 
(r* • "v) fr x "vl 2 


i/2 - 


= T 


£ v 

m X 


r 2 r 


m X 


(67) 


Equations (57), (59), and (67), along with equations (1) and (5) now define an 
optimal three-dimensional trajectory. All of the advantages that were iterated in the 
planar case are still valid; i.e., there are but two adjoint equations which need be 
integrated. The nonplanar component of the steering vector, given by equation (59), is 
carried algebraically, so no further integration is involved in the three-dimensionalization 
of the problem. 


The previously stated disadvantage of r • v appearing in the denominator of 
equation (67) has returned. The obvious need is for a regularization akin to equation 
(18). But in this case such a regularization has not been found. The primary reason for 
this is, simply, the algebraic complexities which result when one attempts to deal with 
the three-dimensional equations. 


A fairly natural approach to avoid the singularity without the difficulties of 
excessive algebraic manipulations is to modify the fundamental representation of 

the X vector, i.e., equation (50). X could be assumed to be written with a vector basis 
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formed by the sets [v, r x (r x v), (r x v)] , [r, v x (r x v), (r x v)] , [r, r x (r x v), (r x v)] , or 
[T, x (? x v), (r x "v)] just as equation (50) assumed (r, v, r x v) as a spanning set. In each 
case the reciprocal r • v occurs, however, so this approach was not successful. 

It is hoped that the difficulty of the singularity will be removed in future research. 

CONCLUSIONS 

It has been known for many years that integrals of the minimum mass 
expenditure optimal trajectory about a single gravitating center exist. These integrals have 
existed as a separate part of the literature and no direct use has been made of them. 

In this paper it has been demonstrated that the use of these integrals yields a 
reduced set of adjoint equations which is to be integrated. This reduced the time required 
for two-point boundary value j-.oolem iterations. 

In the planar case, a regularization has been found which improves numerical 
properties of the equations after the inclusion of known integrals. Such a regularization 
has not been found in the three-dimensional case. 

George C. Marshall Space Flight Center 

National Aeronautics and Space Administration 

Marshall Space Flight Center, Alabama, November 1973 
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